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Abstract 

We consider charged Higgs and W gauge boson contributions to the quasielastic 
scattering + n +p and + p — > r+ + n . These effects modify the standard 

model cross section for these processes and thus impact the extraction of the neutrino 
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standard model cross section, can be significant and can depend on the energy of 
the neutrino. 
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1 Introduction 



We now know that neutrinos have masses and that there is a leptonic mixing matrix 
just as there is a quark mixing matrix. This fact has been firmly estabhshed through 
a variety of solar, atmospheric, and terrestrial neutrino oscillation experiments. The 
next phase in the neutrino physics program is the precision measurement of the 
mixing angles, finding evidence for CP violation and measuring the absolute masses 
to resolve the mass hierarchy problem. 

The existence of neutrino masses and mixing requires physics beyond the stan- 
dard model (SM). Hence it is not unexpected that neutrinos could have non-standard 
interactions (NSI). The effects of NSI have been widely considered in neutrino phe- 
nomenology [Il01SlElEl0SEirara[ia[ia[iatt5in6l[I71[l8l[T9lEOlEIl[2a 
[23l|2aE5l|26l[271|28l|29l[30lEIl[32]. 

It has been established that NSI cannot be an explanation for the standard 
oscillation phenomena, but it may be present as a subleading effect. Many NSI 
involve flavor changing neutral current or charged current lepton flavor violating 
processes. In this paper we consider charged current interactions involving a charged 
Higgs and a W gauge boson in the quasielastic scattering processes Pr + n ^ t~ + p 
and z/t- + p — )■ -f- n. In neutrino experiments, to measure the mixing angle the 
neutrino-nucleus interaction is assumed to be SM-like. If there is a charged Higgs 
or a W contribution to this interaction, then there will be an error in the extracted 
mixing angle. We will calculate the error in the extracted mixing angle. 

The reaction i/^ — )• r~ +p is relevant for experiments like Super-Kamiokande 
(Super-K) El] and OPERA [35] that seek to measure — )• z/^ oscillation by the 
observation of the r lepton . The above interaction is also important for the DONuT 
experiment [SB] which measured the charged- current (CC) interaction cross section 
of the tau neutrino. The DONuT central-value results for a Vr scattering cross sec- 
tion show deviation from the standard model predictions by about 40% but with 
large experimental errors; thus, the measurements are consistent with the standard 
model. The new physics (NP) effects calculated in this paper modify the SM cross 
sections by less than 10% and are therefore consistent with the DONuT measure- 
ments. There have been recent measurements of the appearance of atmospheric tau 
neutrinos by Super-K [33] and by the OPERA Collaboration [35j . 

The reactor neutrino experiments such as Double Chooz [37], Daya Bay [3S] . 
and RENO [39] measure the mixing angle ^13 from the survival probability of an 
electron antineutrino, P( Ve t'e)- If these experiments are designed to observe 
Ve oscillation, the scattering process z/^ -|- p — )■ -|- n will be important. The 

two processes z/^ + ^ — ^ + p and z/^ + p — t- r+ + n involve the same new physics 
(NP) operators. 

Generally, neutrino scattering contains contributions from various processes such 
as quasielastic scattering (QE) , resonance scattering (RES), and deep inelastic 
scattering (DIS). Just above the threshold energy for r production, which is 3.45 
GeV [331 El], the quasielastic interaction dominates in scattering [IHl E]. At 



1 



higher scattering energies other processes have to be included. For instance, the 
DIS is expected to be dominant above around 10 GeV [H], and so Vr scattering at 
the OPERA experiment, running at the average neutrino energy Ey = 17 GeV [35], 
will be dominated by DIS. In this paper we will study the effects of NSI only in QE 
scattering, so we will limit ourselves to energies where QE is dominant. A full NSI 
analysis including all processes in Vr scattering will be discussed in Ref. [42j . 

There are several reasons to consider NSI involving the (z/^, r) sector. First, 
the third generation may be more sensitive to new physics effects because of their 
larger masses. As an example, in certain versions of the two Higgs doublet models 
(2HDM) the couplings of the new Higgs bosons are proportional to the masses, 
and so new physics effects are more pronounced for the third generation. Second, 
the constraints on NP involving the third generation leptons are somewhat weaker, 
allowing for larger new physics effects. Interestingly, the branching ratio of B decays 
to T final states shows some tension with the SM predictions ^31 HI] and this could 
indicate NP, possibly in the scalar or gauge boson sector [15]. Some examples of 
work that deals with NSI at the detector, though not necessarily involving the third 
family leptons, can be found in Refs. [151 HTl US]. 

If there is NP involving the third generation leptons, one can search for it in B 
decays such as B ^ tu^-, B — )■ D^*^tUt- [49j, b — )■ sr+r~ etc. In general, the NP 
interaction in B decays may not be related to the one in Ur + n ^ t~ + p and 
z/^ + J9 — > r+ +n, and so these scattering processes probe different NP. The same NP 
in z^T- + n — 7- + p and z/,- + p — )■ r"*" + can be probed in r decays [50], and we will 
consider the constraint on NP from this decay. However, in general, the scattering 
and the decay processes probe NP in different energy regions. 

The form of NP in z/,- + n — t- + p involves the operator O^p = uTidfV jV^-, 
where Pjj are some Dirac structures. The process v^-^p — > r++n gets a contribution 
from Ojvp- ^"^^ assume CP conserving NP in this paper, and so the coefficients 
of the NP operators are real. The same NP operator can also contribute to hadronic 
tau decays r~ vr^z/^ and r~ — )■ p~z/^, and the measured branching ratio of these 
decays can be used to constrain the couplings in the operator O^p. The ratio of the 
charged Higgs contribution to the SM in z/^ + n — )■ + p and z/,- + p — t- r+ + n is 
roughly (mAr/m^) larger compared to the same ratio in r~ — )■ tt^Ut, where mN,TT are 
the nucleon and pion masses. Hence, significant charged Higgs effects are possible in 
Ur + n^ +p and Ur+p—^T^+n even after imposing constraints from r decays. 
We note that new interactions in the up and down quark sectors can be constrained 
if one assumes CKM unitarity. However, we do not consider this constraint as the 
NP in Onp involves contributions from both the quark and the lepton sectors. 

As noted above, at the quark level NSI in z/^ + n — )■ r~ + p and z/^ + p — )■ r+ + 
involve the u and the d quarks. Often in the analysis of NSI, hadronization effects 
of the quarks via form factors are not included. As we show in our calculation the 
form factors play an important role in the energy dependence of the NP effects. In 
an accurate analysis one should also include nuclear physics effects which take into 
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account the fact that the neutron and the proton are not free but bound in the 
nucleus. There is a certain amount of model dependence in this part of the analysis 
[5T] . and therefore we will not include nuclear effects in our calculation. Such effects 
can be easily incorporated once the free scattering cross sections are known. 

The paper is organized in the following way. In the next section, we present 
a model-independent analysis of NP effects. In the following three sections, we 
consider the neutrino and antineutrino quasielastic scattering in the SM, in a model 
with charged Higgs and in a model with an extra W gauge boson. In the last 
section, we present our conclusions. 



2 Model-independent analysis of new physics 

The process 1/^+12 — ?► +p will impact the measurement of the oscillation probabil- 
ity for the — z/,- transition and hence the extraction of the mixing angle 623- The 
measurement of the atmospheric mixing angle 623 relies on the following relationship 

N{Ur) = P{Uf, -> Ur) X $(zy^J X (TsM(z^r) , (1) 

where N{vt) is the number of observed events, ^{Vfj) is the flux of muon neutrinos 
at the detector, a^^{vr) is the total cross section of tau neutrino interactions with 
nucleons in the SM at the detector, and P^i^^ Vt) is the probability for the 
flavor transition — )■ z/^. This probability is a function of {E, L, Amfj, 6ij) with 
i,j = 1, 2, 3, where Arnfj is the squared-mass difference, 6ij is the mixing angle, E is 
the energy of neutrinos, and L is the distance traveled by neutrinos. The dominant 
term of the probability is 

P(z/^ ^ Ur) ^ sin^ 2^23 cos^ sm^{Aml^L/AE). (2) 

In the presence of NP, Eq. [1] is modified as 

N{ur) = P{u^ ^ Ur) X X atoti'^r), (3) 

with atoti^r) = o'smI'^t) + crNp('^r), where o"np(^'t) refers to the additional terms 
of the SM contribution towards the total cross section. Hence, o"np(^'t) includes 
contributions from both the SM and NP interference amplitudes, and the pure NP 
amplitude. From Eqs. ([Il[3]), assuming 6*13 to be smallp 

sin^ 2(^23) = sin2 2(^23)5A/-^, (4) 

1 + r23 

where ^23 = {(^23) sm + ^23 is the actual atmospheric mixing angle, whereas {023)3 m 
is the extracted mixing angle assuming the SM Uj. scattering cross section. Here 

'^The presence of NP impacts the extraction of the combination sin^ 26*23 cos^ 6*13 . The NP 
changes the extracted value of 623 as well as ^13. But we fix the value of ^13 as an input at this 
point. 
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f2z = o-Mpiyr) I cysMiyT) is the ratio between the NP contribution, including the 
pure NP and interference terms, to the SM cross section. This means that 
can be positive or negative. Assuming neghgible new physics effects in the fi — N 
interaction, the actual mixing angle 623 is the same as the mixing angle extracted 
from the survival probability P(z/^ — )■ z/^) measurement. We will take the best-fit 
value for the mixing angle to be given by 623 = 42.8° [53]. In other words, the 
presence of new physics in a z/T--nucleon scattering will result in the mixing angle, 
extracted from a z/^ appearance experiment, being different than the mixing angle 
from survival probability measurements. The relationship between the ratio r23 
and can be expressed in a model-independent form as 

^ r sin2(6'23)gM 1^ 
'''' ~ [sin2((^23)5Af + 523)J ■ ^ ' 

The reactor neutrino experiments can determine the mixing angle 6*13 from the 
oscillation probability, P(z/e — )■ z/g). The probability of the tau antineutrino appear- 
ance Ue can be used to extract ^13. In this case the effect of NP contributions 
to the process z/,- + p — t- r+ + n is pertinent. The relationship used in measuring 6*13 
will be given as 

N{V^) = P{Ve -> ^r) X $(z>e) X (Ttotl^^r) , (6) 

where [53 ESI |56] 

P{ve Vr) - sin^ 2^13 cos^ 623 sin^ {/^ml^L / AE) . (7) 

Thus the relationship between the ratio of the NP contribution to the SM cross 
section ri3 = crNpipr) / (^sMipT)) which can be positive or negative, and 813 can be 
obtained in a model-independent form as 



ri3 



sin2( 



n3)SM 



1 



Lsin2((6'i3)sM + ^13) 

In Fig. [1] we show the correlation between r23(i3)% and 523(13) [Deg]. One can 
see that ^23 ~ —5° requires r23 ~ 5%. But ^13 ~ —1° requires ri3 ~ 25%. In the 
following sections, we consider specific models of NP to calculate r23 and ri^. We 
will consider a model with a charged Higgs and a W model with both left- and 
right-handed couplings. 



3 Quasielastic neutrino scattering off free nucleon 
- the SM 

In this section we consider the SM contribution to z/^ + n — )■ + p and Vr + P ^ 
r+ + n. We first summarize the SM results for the quasielastic scattering of a 
neutrino on a free neutron target, 

vi{k)+n{p)^l~{k')+p{p'), (9) 
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Figure 1: Correlation plot for = ctj^p^Ut-) / asM{.^T)% versus 52'i[Deg\, and ri3 
(Jnp{,i^t)/(Jsm{,^t)% versus SislDeg]. 



where k, k', p, and p' denote the four- moment a and / indicates the lepton e, /x, or r. 
The spin-averaged matrix element squared for the above reaction is a convolution 
of spin- averaged leptonic and hadronic tensors L'^" and if'^^: 



(10) 



The leptonic tensor calculation is straightforward, but the hadronic tensor involves 
nonperturbative effects. In order to calculate the hadronic tensor, we define the 
charged hadronic current for this process: 



{p{p')\J^Hp)) = v., {p{p')\{V, - A,)\n{p)) 
= VudPip')'i^f,n{p). 



(11) 



The expressions for the matrix elements of the vector and axial-vector currents are 
summarized in terms of six form factors in the Appendix. Due to time reversal 
invariance, the form factors are real functions of t = q^. When invariance under 
charge conjugation holds, two form factors vanish {Fs = 0, Ft = 0) [57j. The 
matrix element, then, can be written as 

M = ^ ' Ui{k')j^{l--f5)u,,{k) 



V2 

un'{p') 



Fr{t)^,+F^{tY''''''' 



2^ +i^^W7.75 + Fp(t)75^ 



M^(p),(12) 



where N and A^' are the initial and final nucleons, while / and vi are the final charged 
lepton and the initial neutrino. In our case N = n, N' = p, I = t, and z/; = i/^-. 

After evaluating |A^P, one can obtain the SM differential cross section for the 
reaction in Eq. ([9]) |57] . 



dt 



M^Gl cos2 



Asm + B 



u] 



SM- 



M2 



SM- 



[ S - U) 

M4 



(13) 
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where Gp = 1.116637 x 10~^ GeV~^ is the Fermi couphng constant, cos^c = 0.9746 
is the cosine of the Cabibbo angle, Mw is the W boson mass, and is the incident 
neutrino energy. M = [Mp + M„)/2 ^ 938.9 MeV is the nucleon mass, and we 
neglect the proton-neutron mass difference. The expressions for the coefficients 
fsM if = A, B, C) are summarized in the Appendix. The Mandelstam variables are 
defined by s = (fc + p)^, t = = {k — k'Y, and u = {k — p')^. The expressions for 
these variables in terms of E^, and the lepton energy Ei are given in the Appendix. 
The quasielastic scattering of an antineutrino on a free nucleon is given by 

ui{k)+p{p) ^l+{k')+n{p'). (14) 

The charged hadronic current becomes |10l |58] 

{n{p')\j;\p{p)) = {p{p)\J^\n{p'))^ 

= Vudn{p)tf,p{p), (15) 

where 

f/.(p,p') = 7ori(p',p)7o- (16) 

The relationship between the differential cross sections of z/^ + n — )■ r~ + p and 
i/^ + p — )■ r"*" + n is 



-^is,t,u) = ^^{u,t,s). (17) 

Thus, the matrix element is given by Eq. [T21 and the differential cross section, 
similarly to Eq. [131 is given by 

dasMii^i) M^Glcos^e,,[ , „ (s-u) ^ (s-u^^ 



dt StxEI 



^SM — BsM — T7o 1" CsM 



M2 



The negative sign of Bsm leads to a relatively smaller cross section for the antineu- 
trino scattering. 



4 Quasielastic neutrino scattering off free nucleon 
— Charged Higgs Effect 

We consider here the charged Higgs contribution to z/,- + n — t- + p and z/^ + p — > 
+ n. Charged Higgs particles appear in multi-Higgs models. In the SM the Higgs 
couples to the fermion masses, but in a general multi-Higgs model the charged Higgs 
may not couple to the mass. What is true in most models is that the coupling of the 
charged Higgs to the leptons is no longer universal. Hence, the extraction of and 
6'i3 from — )■ and z/g — )■ z/g survival probabilities, respectively, will be different 
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from z/^ — )• Vt and z/g — )■ i^,- probabilities, respectively, in the presence of a charged 
Higgs effect. 

The most general coupling of the charged Higgs is 



2^2 



where Ui and dj refer to up and down type quarks, and z/j and Ij refer to neutrinos 
and charged leptons. The other parameters are as follows: g = e/ sm9w is the SM 
weak coupling constant, Vu^dj is the CKM matrix element, and gs,p are the scalar 
and pseudoscalar couplings of the charged Higgs to fermions. Here, in this work, we 
assume the couplings gs^p are real. 

We will choose the couplings gs,p, relevant for z/^+n — i- r^+p and i^r+P — >■ T++n, 
to be given by the two Higgs doublet model of type H (2HDM II). In the 2HDM II 
these couphngs are related to couplings in other sectors and so can be constrained by 
measurements in these other sectors. However, in our analysis, to keep things general 
we will not assume any relation between the couplings gs,p and the couplings in other 
sectors, thereby avoiding constraints from other sectors. To constrain the couplings 
gs,p we will only consider processes that are generated by Onp = uTidfTjUr. In 
the 2HDM II, constrains on the model parameters come from various sectors |60j . 
These constraints turn out to be similar but slightly stronger than the ones obtained 
in our analysis. 

The coupling of charged Higgs boson {H^) interactions to a SM fermion in the 
2HDM II is [61] 

C = ^ cot /3 UiVijPL,Rdj+md^ tBxi(5 UiVijPn^Ldj+m^ tan/3 UiPR^Jj 

(20) 

where Pl,r = (1 -F 7^)/2, and tan (5 is the ratio between the two vacuum expectation 
values (vev's) of the two Higgs doublets. Comparing Eq. f|T9|) and Eq. fl20|) . one can 
obtain 



Uidj 

gs 



rudi tan /3 + m„- cot /3 



Ml 



w 



. _ f rrid^ tan /3 - cot (5 

gp — 



M, 



w 



,,./, ,,,/. m^.tan/? , , 

gt = gp'= ■ (21) 



w 



Constraints on the size of the operator O^p = uVidfVjVj. can be obtained from 
the branching ratio of the decay — )■ vr^z/^. In the presence of a charged Higgs, 
the branching ratio for this process is 
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where the charged Higgs contribution is 



' H 



m^tan^ /3\ m: 



rriu + rud 



m 



(23) 



H 



The SM branching ratio is related to the tau lepton width (r^-) and the decay 
TsXeiV^^^ _ ) as Br^^' _ = T^^f _ with 



16^ 



V mi J 



(24) 



Here 5^/^ = 1.0016 ±0.0014 [62] is the radiative correction. Furtlier, the SM branch- 
ing ratio can also be expressed as [63] 



Brl^^^^-^^ = 0.6075r(r- ^ u^e-^^ 



10.82 ±0.02% 



(25) 



while the measured i?r(r — )■ tt z/r)ea;p = (10.91±0.07)% [61]. In Fig.[2]we show the 
constraints on m/^ — tan /3 from — t- tt^z/^. From Eq. |20]we can construct the NSI 
parameters defined in Ref [H] as €rr^^^ = and = . We find 



that the constraints on the effective operator considered in this work are consistent 
with the one in Ref. [IS]- Finally, we note that r has a significant branching ratio 
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Figure 2: Constraint by Br{T — )■ tt z/^) at 95 % CL. The colored region is allowed. 



to —7- p^i'r [64] • However, a charged Higgs cannot contribute to this decay, and 
hence there is no constraint on the charged Higgs couplings from this decay [50] . 

Keeping in mind the constraints from Fig. ([2]), we calculate the charged Higgs 
contribution to z/^ ± n — )■ r~ ± p. The modified differential cross section for the 
reaction in Eq. (j9]) is 



da 



SM+H 

~dt 



M^Gl cos2 
87rE2 



Ah + B 



H- 



M2 



c. 



SM- 



[s — u 
M4 



(26) 
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where xh 



m 



w 



/Ml, Ah = Asm + 2xHi?e(A^) + x^A^, and Bh = Bsm + 



2xhRg{B\j). Superscripts I and P denote the SM-Higgs interference and pure 
Higgs contributions, respectively. The expressions for the quantities A^j^ and B\j 
are given in the Appendix. The terms A^ and B\j are proportional to the tiny 
neutrino mass, and we will ignore them in our calculation. Note that this happens 
because we have chosen the couplings to be given by the 2HDM II. With general 
couplings of the charged Higgs, these interference terms will be present. The charged 
Higgs contribution relative to the SM r^^ = asM{y\ proportional to t because of 
the dominant term XtC^,, where Xt = t/4M'^ (see the Appendix for more details). 
Consequently, is proportional to the incident neutrino energy (see Fig. (j3])). The 
deviation is negative, as there is no interference with the SM; hence, the cross 
section for I'r + n ^ t~ + p is always larger than the SM cross section. This means 
that, if the actual 623 is close to maximal, then experiments should measure 623 
larger than the maximal value in the presence of a charged Higgs contribution. 

The differential cross section for the interaction z/^ + j9 — )• r+ + n has the same 
form as Eq. [2n]in the limit of a massless neutrino. The hadronic current in this case 



is the complex conjugate of the one in the Appendix. The ratio — ^^^^(p y 
well as the deviation ^13, is shown in Fig. HI As 613 is a small angle, large tan /3 and 
small charged Higgs mass are preferred to produce an observable deviation 613. For 
instance, we find S13 ^ 1° and r]f ^ 30% at = 8 GeV, Mh = 200 GeV, and 
tan/3 = 100. 



as 
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Figure 3: Variation of r]f% with E,y and variation of ^23 with Mh and E,y. The 
green line corresponds to the SM prediction. The black (dotdashed), pink (dashed), 
and blue (solid) lines correspond to tan/3 = 40, 50, 60. The right figure is evaluated 
ai Ely = 5 GeV, while the left figures are evaluated at Mh = 200 GeV. Here, we use 
the best-fit value ^23 = 42.8° |53]. 
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Figure 4: Variation of r)^% with E^, and the variation of (5i3 with Mh and E^. The 
green hne corresponds to the SM prediction. The black (dotdashed), pink (dashed), 
and blue (solid) lines correspond to tan/3 = 80, 90, 100. The right figure is evaluated 
&i Ey = GeV, while the left figures are evaluated at Mh = 200 GeV. Here, we use 
the inverted hierarchy value 6*13 = 9.1° 

5 Quasielastic neutrino scattering off free nucleon 
- W model 

Many extensions of the SM contain a W gauge boson. We next consider modification 
to Vr + n ^ + p and z/^ + p — t- r+ + n in models with a W . There are limits 
on the W mass from direct searches to final states involving an electron and muon 
assuming SM couplings for the W [61]. These limits generally do not apply to the 
W coupling to and r which is relevant for our calculation. 

The lowest dimension effective Lagrangian of W interactions to the SM fermions 
has the form 

C = ^VfjfY{9pPL + 9{^^PR)fW'^+ h.c, (27) 

where / and / refer to the fermions and g)^ )j. are the left- and the right-handed 
couplings of the W. For a SM-like W boson, g{^^ = 1 and g^^^ = 0. We will 
assume g{^ to be real. Constraints on the couplings in Eq. ( p7|) come from the 
hadronic r decays. We will consider constraints from the decays r~ tt^Uj. and 

The branching ratio for — )■ n^Ur is 

5rf"i^:. = Br'JU.-..il + rw')\ (28) 
where the W contribution is 

rw' = xw'g7{gt-gt), (29) 

and xw' = ^w/M'^,. The branching ratio for the — )• p^Ur process is 

Bri^::-:^ = 5ev..(l+-^')^ (30) 
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with the W contribution 

r^. = xwgTigr + g'-n')- (31) 

The SM branchine ratio is related to the decay rate as Br^!^\ _ = Y^^'\ _ /Tr 
with 



G 



(32) 



where L = 223 MeV 



. Further, the SM branching ratio can also be expressed 



as 



Br, 



SM 



1.235r( 



r — 7- z/^e z/g 



21.92 ±0.05%. 



(33) 



The measured branching ratio is Br{T^ — )• p^Vr)exp = (23.1 ± 0.98)% 

Figures |5] and [6] show the allowed regions for the W couplings. The couplings are 
uniformly varied in the range [—2,2] and constrained by the measured — )■ n^Vr 
and T~ — )■ p~Vr branching ratios with la errors. From Eqs. fj29| |3T]) . the case with 
a pure left-handed W coupling is allowed, as shown in Figs. ([5l[6]). The constraints 
on the effective operator are consistent with the one in Ref. [H]. From Eq. [271 the 



NSI parameter e 



ud{L,R) 



defined in Ref. 



ud(L,R) 

IS given as Ett 



„TV ^ud ( Mw\2 

9l 9(^l,r)\m^,) ■ 
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Figure 5: The constraints on the W couplings without right-handed coupling at 
My/' = 500 — 1000 GeV. The constraints are from r~ — )■ tt'i/t- and r~ — > p'^^ 
branching ratios. The errors in the branching ratios are varied within la. The 
colored regions are allowed. 

In the presence of the W gauge boson, we can obtain the modified differential 
cross section for the reaction z/^ + — t- r~ + p as 



Jt 



M^Gl cos2 



u 



M2 



M4 



(34) 
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Figure 6: The constraints on the W couphngs with both left- and right-handed 
couphngs at Mw' = 500 — 1000 GeV. The constraints are from — t- tv^Ut and 
— )■ p'l'r branching ratios. The errors in the branching ratios are varied within 
la. The colored regions are allowed. 



where the coefficients A', B', C include both the SM and W contributions. The 
expressions for these coefficients are given in the Appendix. 

For a SM-like W boson, with right-handed couplings ignored, the structure of 
the differential cross section is similar to the one in the SM case. Hence, the W 
contribution relative to the SM r??r, = '^'^'^'^^\ does not depend on the incident 

neutrino energy E^. We find r^, ~ 5% at Mw' = 500 GeV from the hadronic tau 
decay constraints in Fig. ([5]). The variation of ^23 with the W mass is shown in 
Fig. ([7]). In this case, is always negative and can reach up to —5° at Miy, = 500 
GeV. Note that ^23 does not dependent on E^, either. 




700 750 800 850 900 950 1000 '4 5 6 7 8 

Mw [GeV] Ey [GeV] 



Figure 7: The left (right) panel illustrates the deviation ^23 with the W mass 
{E„) when only left-handed W couplings are present. The lines show predictions 
for some representative values of the W couplings {gY^idlf") taken from Fig. (|5]). 
The green line corresponds to the SM prediction. The blue (solid, lower) line in the 
left figure corresponds to (0.69, 0.89) at = 5 GeV, and the blue (solid, lower) 
line in the right figure corresponds to (1.42, 0.22) at M\y' = 500 GeV. Here, we use 
the best-fit value ^23 = 42.8° [53j. 
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Next, we consider the right-handed couphngs also. The variation of r^,% with 
Mw' in this case is shown in Fig. (|8]). The r^,% values are mostly positive which, 
in turn, leads to ^23 being mostly negative. We find that r^,% depends slightly on 
the neutrino energy. The variation of ^23 with the W mass and Ei, are shown in 
Fig. ©. 

The W contribution to the interaction z/^ + p — > r+ + n leads to the following 
differential cross section: 



dt 



87rE2 



A' -B 



M2 



2-, 



M4 



(35) 



The differential cross section of the antineutrino scattering is relatively smaller than 
the corresponding one for the neutrino scattering because of the negative sign of 



w 



the B coefficient. Thus, the value of the ratio 
corresponding ratio. 



,23 
W 



is smaller than the 



Because of the smallness of 613 and rj^,, the NP effect on 



the extraction of 6*13 is small. Achieving large rl^, within the constraints given in 
Fig. [6] is difficult in this model. This means the effect of the NP contribution in 613 
is very small and we do not plot the results of this calculations. 



^2 



O F 

500 600 



700 800 900 1000 
Mw [GeV] 



5 6 
Ey [GeV] 



Figure 8: The left (right) panel illustrates the variation of r^,% in Ur + n ^ t~ +p 
scattering with the W mass {E,y) when both left- and right-handed W couplings 
are present. The lines show predictions for some representative values of the W 
couplings {qY"^ 1 9^ 1 dlii) taken from Fig. The green line corresponds to the SM 
prediction. The blue (solid, upper) line in the left figure corresponds to (-0.94 , 
-1.13 , -0.85) aX = GeV, and the blue (solid, upper) line in the right figure 
corresponds to (1.23 , 0.84 , 0.61) at Mw' = 500 GeV. Here, we use the best-fit 
value ^23 = 42.8° |53]. 
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Figure 9: The left (right) panel illustrates the deviation ^23 with the W mass (E,^) 
when both the left- and right-handed W couplings are present. The lines show 
predictions for some representative values of the W couplings [gY'' , , g^) taken 
from Fig. dH]). The green line corresponds to the SM prediction. The blue (solid, 
lower) line in the left figure corresponds to (-0.94 , -1.13 , -0.85) aX = GeV, 
and the blue (solid, lower) line in the right figure corresponds to (1.23 , 0.84 , 0.61) 
at Mw> = 500 GeV. Here, we use the best-fit value 623 = 42.8° [M]. 



6 Conclusion 

In this paper we calculated the effect of a charged Higgs and a W contribution to 
Ur + n ^ T~ + p and z/,- + p — )■ r"*" + scattering. We constrained the parameters 
of both the models from — )■ tt^z/^ and r~ — )■ p~z/^ decays. Corrections to the 
SM contribution to Ur + n ^ -\- p and + p ^ t"^ + n impact the extraction 
of the neutrino atmospheric mixing angles 623 and 613, respectively. We found that 
the charged Higgs model can produce significant corrections to ^23, 13 that measure 
the deviation of the actual 6'23,i3 from the (6'23,i3)5a/ angles which are extracted 
assuming the SM v^-lvj. scattering cross sections. The W model effect generates a 
large deviation ^23 but negligibly small 813. As 613 is smaller than ^23 larger NP in 
z/r + P — ^ + ^ is required to produce effects in 813 similar in size to 523. 

When a charged Higgs is involved, ^23,13 are negative. This is because there 
is no interference of the charged Higgs contribution with the SM contribution, for 
massless neutrinos, and so the cross sections for z/^ + n — )■ t~ +p and z^r+P — ^ +n 
are always larger than the SM cross sections. This means that experiments should 
measure ^23,13 larger than the present values in the presence of a charged Higgs 
contribution. In this case we also found that ^23,13 increase in magnitude with the 
neutrino energy. Hence, a possible sign of the charged Higgs effect would be a 
measurement of ^23,13 that shows an increase with increasing neutrino energy. 

For the W model we calculated a significant contribution to ^23 which can be 
both positive and negative, but is mostly negative. The deviation ^23 was found 
to be independent of the neutrino energy for a left-handed W but neutrino energy 
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dependent when both left- and right-handed W chirahties were present. A neghgibly 
small deviation, 613, was found in the W model because of the small value of 613. 

We have presented in this paper a first estimation of the charged Higgs and W 
effects in the extraction of 6^23 and 6*13. We hope more detailed calculations including 
nuclear as well as detector effects, will be done to find out whether these new physics 
effects can be observed at present Vt/^t appearance experiments and/or to motivate 
new experiments that can detect these effects. 
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Appendix 

Hadronic form factors 

The expressions for the vector and axial- vector hadronic currents in Eq. SH are 

{p{p')\yM{p)) = Up{p') 

-{p{p')\Af,\n{p)) = Up{p') 

Here q = p' —p and the form factors Fi are functions of t = g^. The parametrizations 
of the axial- vector and pseudoscalar form factors are [57] 

FA{t) = ^a(0)(i-^)"', 

Pp(t) = 37 

where Fa{0) = —1.2695 is the axial coupling |i64], ""^tt is the charged pion mass, 
and Ma = 1.35 GeV is the axial- vector mass [67]. The expression for Fp{t) can be 
shown to be true at low energy, where the predictions of chiral perturbation theory 
are valid [68]. We have assumed the relation to hold at high t also. Note that Fa{0) 
is sometimes replaced by FA{t), which gives similar results for Fp{t) at low t but 
very different results at high t. 

The Dirac and Pauli form factors FY2 ^^^^ 

TpVu\ G E{t) - XtG M{t) T?V(+\ GM{t)-GE{t) , . 

j'lW = Yzr^^ ' ^2 w = — ^3^7 — ' 1^8) 



M 



l5Un{p)- (36) 
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where Xt = t/AM"^ and 

Kjm — ^ M ~ — ^ E ~ ^ E ■ V'^) 

Here G^" and are the electric and magnetic form factors of the proton and 
neutron, respectively. The simplest parametrizations of these form factors are given 
by the dipole approximation 

G\ ^ Gd, Gl ^ 0, Gl, ^ fipGn, G^/ ^ f^nG^ , (40) 

where Go = (1 — t/My)'^, My = 0.843 GeV is the vector mass, and fJ^pif^n) is the 
anomalous magnetic moment of the proton (neutron) [69j. 

In the presence of the charged Higgs, applying the equation of motion to the 
hadronic matrix elements for the scalar and pseudoscalar currents for the process 
Ut- + n ^ T" + p gives 

{pip'Mgp - gf^i,)d\nip)) = Pip,) {gT'Gs + 9p''Gp^,) nip,), (41) 



or 



{pip)\ud\nip)) = pip4.)Gsnip2), 
-{Pip'))\ui5d\nip)) = p(p4)Gp75n(p2) , (42) 



where 



Gsit) = r^FYit), with = ^\ ~ 0(1), 

TUg 

with rhg = irriu + mrf)/2. In the W model, the current has both V ± A structures. 
One has to calculate the matrix element, 

(p(p')l j;i^(p)) = {Pip')\u {gfi.il - 75) + 9R%il + 75)) d\nip)). (44) 
Kinematic details 

The Mandelstam variables in terms of and the lepton energy Ei are 

s = + 2ME^, t = 2MiEi- E^), 
s-u = AME^ + t-ml (45) 

Then t and Ei lie in the intervals 

- 2El^ iEr + pD <t<m]~ 2ET iEr - pD , (46) 
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2M _ i _ ^ ^ 2M > V ; 

where the energy and momentum of the lepton and the neutrino in the center of 
mass (cm) system are 

i s - M^ ) 
2v^ 



The threshold neutrino energy to create the charged lepton partner is given by 

where m;, Mp, M„ are the masses of the charged lepton, proton, and neutron, 
respectively. In our case, the threshold energy of the tau neutrino is i?*^ = 3.45 
GeV. 

The differential cross section in the laboratory frame is given by 

d(ytoM) ^ \M\'' 

dt U-kEW' ^ ' 



The expressions for the coefficients fsM (/ = -B, C) in the SM differential 
cross section [see Eq. f|T3|) ] are 

Asm = 4(xt - xi) [{E^)\l + xi + Xt) + (F^)'(-l +xi + Xt) + {E^)\xi + x\ + Xt) 
+AFlxiXt + 2F^F^{xi + 2xt) + AFaEpXi^ , 

BsM = AxtFAiF^ + F^) , 

{FYf + Fl-x,{Fyf 

where xi = mf /AM"^. 

The expressions for the quantities A^^f and Bj^ in the differential cross section 
in Eq. (l26|l are 

A'h = 2^i{xt-xi)gp{g's''-g'p)Gp{FA + 2FpXt), 
K = lv^i9s'i9s'-9p)Gs{Fr + F^x,), 

Al = 2(x,-xO(|^7?P+|^7^^f)(|^7]^TG^x, + |^7^tG|(xi-l)). (52) 
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For the 2HDM II model couplings, gs,p are given in Eq. (121 p . Note that the inter- 
ference terms A^^j and Bfj vanish in this model. 

The expressions for the quantities /' (/ = A, B, C) in the differential cross 
section in Eq. are 



A' = A{xt - xi) [(1 + r'^,y(^{F^y{l + XI + Xt) + 2F^F^{xi + 2xt) + {F^)\xi + x\ + Xt)) 

+ (1 + rl^)^{[FAf{-^ + + Xt) + ^FaFpx, + ^Flxixt^ 
B' = 4Re[{l + r'^,){l + r-^,)]xtFA{F^ + F^), 



C 



1 



1 + r'^'YUFrr - ^tiF^f) + (1 + r-w'Yn ■ (53) 



In the absence of W contributions, the /"s reduce to the respective SM results in 
Eq.EI 
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